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We develop a general framework for effective equations of expectation values in 
quantum cosmology and pose for them the quantum Cauchy problem with no-boundary 
and tunneling wavefunctions. We apply this framework in the model with a big negative 
non-minimal coupling, which incorporates a recently proposed low energy (GUT scale) 
mechanism of the quantum origin of the inflationary Universe and study the effects of 
(^ • the quantum inflaton mode. 
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1. Introduction 

Efim Samoilovich Fradkin created a scientific school that can be characterized by a 
peculiar style that combines great diversity of fundamental physical problems and efficient 



X 



methods of their solution. He never worked in the area of quantum cosmology, but the 
generality of his methods, either it is the functional approach to QFT, the Euclidean QFT 
or the most advanced methods of gauge constrained theories, not only apply in this field but 
actually become indispensable. The authors of this paper had a happy opportunity to share 
a small piece of his scientific wisdom and, in this work, would like to show how these three 
methods pioneered by E.S. Fradkin fruitfully combine in quantum cosmological context. 

The power of his approach is based on the fact that the consideration of the problem in 
question begins with the first fundamental principles followed by a set of clear mathematical 
transformations which unambiguously push the solution to its logical extreme. This strat- 
egy straightforwardly leads to the final result instead of unreliable hand waving and (very 
often erroneous) insight. In context of quantum cosmology this strategy implies a concrete 
setting of the Cauchy problem for the initial cosmological state, the unitary evolution from 
which gives rise to the present-day semiclassical Universe with its observable cosmological 
parameters - Hubble constant H, density parameter Q, anisotropy of microwave background, 
etc. Here we show how this strategy can be realized for a particular model of inflationary 
Universe. 
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2. Quantum origin of the Universe as a low-energy 

phenomenon 

The virtue of quantum cosmology is that it can give initial conditions for inflation, 
which in their turn determine main cosmological parameters of the observable Universe, 
including the density parameter Q. The requirement of Q > 1 in closed cosmology gives 
the bound on the e-folding number N > 60 - the logarithmic expansion coefficient for the 
cosmological scale factor a during the inflation stage with a Hubble constant H = a/a. In 
the chaotic inflation model H = H(p) is generated by the inflaton tp and, therefore, all the 
parameters can be found as functions of initial tp. This quantity is subject to the quantum 
distribution p(p) determined by the cosmological wavefunction. If this distribution has a 
sharp probability peak at certain tp = p>i, then this value serves as the initial condition for 
inflation. 

Two known quantum states that lead in the semiclassical regime to the closed inflationary 
Universe are the no-boundary Jl| and tunneling |2| wavefunctions. They both describe quan- 
tum nucleation of the Lorentzian quasi-DeSitter spacetime from the Euclidean hemisphere - 
the gravitational instanton responsible for the classically forbidden state of the gravitational 
field. Tree level wave functions are generally devoid of the observationally justified proba- 
bility peaks. Beyond the tree level the distribution pnb,t(<^) becomes the diagonal element 
of the reduced density matrix obtained by tracing out all degrees of freedom but (p |3|, Q] 

Pnb,t(<^) ~ exp[^/(v?) - r((p)\. (2.1) 

Here the classical action is amended by the Euclidean effective action r(ip) of all quantum 
fields that are integrated out. 

The model of || capable of generating the probability peak contains the graviton-inflaton 
sector with a big negative constant, — £ = |£| ^> 1, of non- minimal curvature coupling, 

L{g^) = 9 1/2 [(^ - \& 2 )R ~ \(Vp) 2 - ~mV - \^) , (2.2) 

and generic GUT sector of Higgs x, vector gauge A^ and spinor fields ip coupled to the 
inflaton via the interaction terms \ x x 2 <p 2 , g\A 2 tp 2 , fyiptpip (non-derivative parts) with the 
coupling constants \ x ,gA,ftp- Its peak-like distribution function can be approximated by 
the gaussian packet 

Pnb,tM - ~n^ ex V f-(v - ^) 2 /2A 2 l , (2.3) 
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where the parameters of the peak - mean value pi = mp(87i\l + 5\/\^\A) 1 ^ 2 and its quantum 
width A = (</?//\/12A)vA/|£| are strongly suppressed by a small ratio vA/|£| known from 
the COBE normalization for AT/T ~ 10~ 5 (because the CMBR anisotropy in this model 

is proportional to this ratio [[7j). Here A = 1/2A f Y, x ^x + ^^a9a ~ 162^/, /$ J ^ s ^ ne 
universal combination of coupling constants above and 

5 = -^\i\m 2 /\m 2 P . (2.4) 



For the no-boundary and tunneling states (±-signs respectively) the peak exists for positive 
A and ±(1 + 5) < 0. The classical equations of motion in the slow roll approximation, 

<p + 3H((p)<p-F(<p) = 0, (2.5) 

H(tp) ~ (A/121^1) 1 / 2 ^, F(ip) ~ — \m 2 P {l + £)y9/487r£ 2 , show that the inflaton decreases from 
its initial value, (p ~ F/H < 0, only for 1 + 5 > 0, that is only for the tunneling state. 
Only in this case the duration of the inflationary epoch is finite with the e-folding number 
|J N — J dtH ~ ASix 2 / A. Comparison with N > 60 yields the bound A ~ 5.5 which 
can be regarded as a selection criterion for particle physics models ||. This conclusion 
remains qualitatively true when taking into account the contribution of the inhomogeneous 
quantum modes to the radiation current of theeffective equations || . This contribution and 
its dynamical effect were obtained in || by the method of the Euclidean effective action, 
however, the quantum fluctuations of the inflaton field itself have not been taken into account. 
For the proponents of the no-boundary quantum states in a long debate on the wave- 



function discord 0, [HJ this situation seems unacceptable. According to this result the 
no-boundary proposal does not generate realistic inflationary scenario, while the tunneling 
state does not satisfy important aesthetic criterion - the universal formulation of the initial 
conditions and dynamics in one concept - spacetime covariant path integral over geometries. 
Thus, one of the motivations of considering the quantum mechanical sector of the inflaton 
mode is the hope that it can handle this difficulty In view of the smallness of A the quantum 
fluctuations A(p ~ A are expected to be negligible, but those of their quantum momenta 
App ~ 1/A blow up for small A. Therefore, apriori, it is hard to predict the overall mag- 
nitude of the quantum rolling force and its sign due to A(p(t). In what follows we carefully 
consider this problem. 

3. Effective equations: setting the problem 

Effective equations for expectation values of quantum fields, g(x) = 0(x), x(x),ip(x), A^(x), g^ u {x) 
in the quantum state |\I/), g(x) = (\E r |g(x)|\E r ), can be obtained by expanding the Heisenberg 
equations of motion for g(x) = g(x) + Ag(x) in powers of Ag(x) and averaging them with 
respect to |\I>): 

g + ■/<») = <>. (3.1) 

Here S[g] is the classical action and J(x) is the radiation current which accumulates all 
quantum corrections starting with the one-loop contribution 
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+ ..., (3.2) 



and G(z,y) = (\&| Ag(z) Ag(y) |\I>) is the Wightman function of quantum disturbances in a 
given quantum state. 

In closed cosmological model, the total metric and the scalar field (playing the role of an 
inflaton) is usually decomposed in the spatially homogeneous background and inhomogeneous 
perturbations ds 2 = — N 2 {t)dt 2 + a 2 {t^ijdx 1 dx^ + h^ v {x)dx^ dx u , f{x) = ip(t)+8ip(x), where 



a(t) is the scale factor, N(t) is the lapse function and 7^ is the spatial metric of the 3-sphere 
of unit radius. Therefore, the full set of fields g(x) consists of the minisuperspace spatially 
homogeneous variables Q(t) and inhomogeneous fields f(x) depending on spatial coordinates 
x*=x, Q(t) = (a(t), <p(t), N{t)), f{x) = {6<p(t,x), V(t,x), x (i,x), ^(*,x), A M (*,x), ...)■ 
From the structure of the initial quantum state, that will be discussed later, it follows that 
only minisuperspace variables have nonvanishing expectation values ( Q(t) ) 7^ 0, ( /(x) ) = 
0. Therefore, the full set of effective equations reduces to the following three equations in 
the minisuperspace sector 

Tq77\ + JQ(t) =0; Jq — Jn, J a , J v , (3.3) 

their quantum radiation currents Jq(£) containing the contribution of quantum fluctuations 
of minisuperspace modes themselves and those of the spatially inhomogeneous fields. 

Splitting of the whole configuration space into minisuperspace and inhomogeneous sectors 
reflects the choice of the collective variables and their quantum states - they turn out to 
be very different for these two sectors of the theory. This results in different calculational 
strategies for the corresponding quantum averages. Inhomogeneous /-modes live in the 
Euclidean DeSitter invariant vacuum, so that their calculation is strongly facilitated by the 



Euclidean effective action method |l .1J. For the quantum mechanical <£>-mode the situation 
is different. 

In the tree-level approximation the wavefunction exp[=p/(<^)/2] of (p does not have good 
probability peaks and is even unnormalizable. Therefore, the tree-level quantum average 
(AQAQ) rie is badly defined. Beyond the tree-level approximation the situation can be 
improved, because the quantum average should now be defined with the aid of the reduced 
density matrix (AQAQ) = tr AQAQ p , 

p = p(cp, <p') = J df*{cp, /)* V, /)> ( 3 - 4 ) 

which originates from tracing the /-variables out and includes loop corrections. As shown 
in H |I2| , the diagonal element of this density matrix - the distribution function of tp, 
p(<p) = p(<p,ip), is given in the approximation of a gaussian integral by the effective action 
algorithm ( |2.1|) which can generate a sharp probability peak ( |2.3|) . With this modification the 
quantum correlators become well defined, being expressed in terms of (AipAip) ~ A 2 < 00. 

4. Quantum Cauchy problem: tree level approxima- 
tion 

In the one-loop approximation the radiation current can be calculated on the classical 
background - the lowest order approximation for the mean field. The initial conditions for 
this solution follow from the no-boundary and tunneling wavef unctions. We show this for the 
model of a minimally coupled inflaton field <ft with a generic potential V((p) (we reserve the 
notation <ft as opposed to the notation (p for the non-minimal inflatonF]). The Wheeler-DeWitt 

^^This framework can be extended to the non-minimal model by reparametrizing the latter to the Einstein 
frame fl3fl, and this will be done below. 



equation for the cosmological wavefunction on 2- dimensional minisuperspace, &(</>, a), has 
two semiclassical solutions - the so-called no-boundary and tunneling wavef unctions. In the 
approximation of the inflationary slow roll (when the derivatives with respect to are much 
smaller than the derivatives with respect to a) they read |§ 



* NB {<t>,a) = C NB {a 2 H 2 {<f)) - 1)-V4 e -W2 cos 
* T (0,a) = C T (a 2 H 2 (<p) - l^V^exp 
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(4.1) 
(4.2) 



and describe the nucleation of the Lorentzian DeSitter spacetime with the Hamilton- Jacobi 
function S((f>, a) = —nm 2 p (a 2 H 2 ((j)) — l) 3 ^ 2 /2H 2 ((p) from the gravitational half-instanton with 
the action J(0)/2 = — 7nn'p/H 2 (<j)). This nucleation takes place at a = 1/H(<f>), H((p) = 
(«;y(0)/3) 1 / 2 being the effective Hubble constant, driving the inflationary dynamics of the 
model, a/ a ~ H(<f)) (k = 87r/m 2 P - the gravitational constant). This domain forms the 
one-dimensional curve in two-dimensional minisuperspace, 

x(<M=a-(3/ K \/(0)) 1 / 2 = O, (4-3) 

which can be identified with the physical subspace, xi'Pi a ) being regarded as the correspond- 
ing gauge condition. If <ft is chosen as a physical coordinate then, according to the formalism 
of JL4], the physical wavefunction Vl/(0) follows from the Dirac wavefunction \&(0, a) by 
the transformation \I/(0) = J 1 / 2 ^^, a)| x (^ a )=o, where J x is the Faddeev-Popov determi- 
nant - the Poisson bracket of the gauge condition with the Hamiltonian constraint. Since 
J 1 / 2 ~ (H 2 a 2 — l) 1 / 4 , it cancels the divergent preexponential factors in (|4.1|)-(|4T^) and results 
in the initial physical wavefunctions ^nb,t{4>) ~ ex P[T-^(0)/2]. Beyond the tree level they 
get replaced by the density matrix ( |3.4| ) which, for the model of the nonminimally coupled 
inflaton, even yields a sharp probability peak and, thus, generates the expectation value of 
the inflaton <fi = (<f)) - the first initial condition of the above type. 

The second initial condition arises from the expectation value of the physical momentum 
conjugated to 0, p^ = (p^). In view of reality of the initial density matrix (|3.4j) this 
expectation value is vanishing, (p<f>) = 0. From the Hamiltonian reduction of the symplectic 
form in the gauge x(0 5 a) = it follows that the physical momentum expresses in terms of 
the original momenta, U a da + li^dcj) = p^dcj), p 4> = Ii <i) - Ii a X<t>/Xa, X<t> = d$X, Xa = d a x- 
Therefore, for p^ = 0, n^ homogeneously expresses in terms of n a and, in view of the 
Hamiltonian constraint this implies that at the initial Cauchy surface n^, = and n a = 0. 
Thus, the full set of initial conditions for the classical background reads as — (</>)> a — 
l/H(<p), = a = O. 

5. Cauchy problem for linearized quantum fluctuations 



In the Hamiltonian reduction to the physical sector for cosmological perturbations |L5 
scalar perturbations of metric and inflaton fields (which only contribute to the homoge- 
neous sector) are defined by the ansatz: ds 2 otal = a 2 (r]) [—(1 + 2A)dr] 2 + (1 — 2ip)jijdx 1 'dx i ] , 
0totai = + S<f). Their quadratic action reduces to the functional of only two invariant (with 
respect to linearized gauge transformations) combinations of canonical coordinates (ip, 8(f) 



and their momenta (IT,/,, ILj^), ^ = ip + 7Y50/0', Uy, = IT,/, — 2a 2 v /7_D50/K0', where prime 
denotes derivatives with respect to the conformal time rj, Ti — a' /a is the conformal time 
Hubble constant and D = 7 U VjVj + 3. With the choice of physical phase space variables, 
q = 2a^ / K(j)' + HD~ 1 U 1 / J /a(j)'y/ ; y, p = —<fi'a y / : yDty/2H + K,(j)'Tl^/4a, the action acquires the 
final Lagrangian form 

S[q}\ 2 =~ J d V ^(-Dq)[- d 2 /d V 2 + cj ) '(l/ ( f ) y + K<f>' 2 /2 + D}q. (5.1) 

The field q is well known from the theory of cosmological perturbations as the Bardeen 
invariant [|3~5j]p|. 

In the Newton gauge, widely used in the theory of cosmological perturbations, these 
perturbations read in terms of the Bardeen invariant q as [I5| A — if), tp — n<f)'q/2a, 



V q)' / a<p' . Similar relations in the minisuperspace gauge ( |4.3| ) can be obtained by 
expressing 5a = —aip in terms of 50, ip = (x<?i/ a Xa)50, and finding the expression for the 
physical momentum p S4> conjugated to 50, Yl^ip'+Yls^cf)' = p S(f ,6 &'+..., p S(/> = n^+IL^/aXa- 
Then, the linearized Hamiltonian constraint gives IT,/, and LT^ as functions of (5(j), ps^), and 
one directly proceeds to the transformation relating (q,p) to (50,^). For 77 — ^ this 
transformation turns out to be singular (see |ll[] for details). This singularity is, however, 



an artifact of the definition of the invariant variables non- analytic at 0' — » 0. To see this, we 
decompose the general solution for q(rj) in the sum of two linearly independent solutions of 
the equation of motion for the action ( p.l[ ), q(rj) = c + q + (j]) + c_q_(r/), one of them having 
a singular behaviour at r\ = 0, q + {rf) = rf/4>' (1 + 0(r] 2 )), q_(r}) = l/4>' (1 — 3r/ 2 /2 + 0(r/ 3 )). 
Then we show that the coefficients c± have a regular solution in terms of the initial conditions 
for physical variables (50(0), p S(f> (0)) Q, c + = -V^(0)/9#V7, c_ = V^<f)(0)/3H 3 / 2 . 
This relation will be used throughout the rest of the paper to express the Heisenberg opera- 
tors A<5 P hys(?7) and AQ(rj) in terms of the Schrodinger ones, 50(0) = <50, p^(0) = d/id(8(f)), 
and then find the quantum averages of their bilinear combinations in the ^-representation 
of the initial density matrix (see eq. Q2.3p with 50 = — (0)). 

6. Non- minimal model 

The action of the model ( |2.2| ) has a generic form 
S[g^if} = Jd 4 xg 1 ' 2 |-V(^) + U(<p)R - ^(^(V^) 2 } , (6.1) 

where the coefficient functions can be read off (|2.2|) . In the presence of spatial densities of 
one-loop radiation currents jjy = Jn/o 3 ^, jtp = J^/Na 3 ^/, the effective equation for the 

2 Note that the operator D is positive definite for the spatially homogeneous mode, D = +3. Thus, this 
is a ghost variable signifying the classical instability. This instability at the linear level is the manifestation 
of inflation which is a huge instability phenomenon incorporating the runaway modes. In contrast with 
the S-matrix theory this instability should not be regarded as an irrecoverable flaw of the theory, because 
we know a nonlinear mechanism that provides a graceful exit from the inflation stage. In particular, no 
special measures like introducing the indefinite metric should be undertaken to eradicate this phenomenon. 
Homogeneous fluctuations of the inflaton field do not have a particle nature and one should not take care of 
guaranteeing the energy positivity of their excitations. Therefore, this mode can and should be quantized in 
the coordinate representation with positive metric in the Hilbcrt space. 

6 



inflaton field reads ITT 



(/3+(3 — U v j^(p-F((p,a, (p) = 0, (6.2) 

% ' a ^ ) = T^4i ln((7 + 3[/ ' )+floop(v ' a '^ d) ' (6 ' 3) 

F loop (^, a, 0, a) = Z7q ^ (W - 2U v3N - ±^fj . (6.4) 

It contains quantum contributions to the friction term and the rolling force, while the first two 
terms in ( |6.3|) represent the classical part. As regards the quantum part, one-loop radiation 
currents split into the contributions of the quantum mechanical sector and the sector of 
spatially inhomogeneous modes, ji_i 00 p = j 9 +j^ ■ The /-part of the current can be generated 
by the effective action, which implies the replacement of the classical coefficient functions 
V(<p), U(if), G(ip), by their effective counterparts and truncation of the (generally infinite) 
series to the first three terms. This truncation is based on two assumptions - the smallness 
of inflaton derivatives due to the slow roll regime and smallness of R/rn?^ - the curvature 
to particle mass squared ratio [§]. Thus, with this approximation, the effective equations 
of motion in our non-minimal model take the form of fl6.2|) with V e s ((p), U e s(p), G e g (ip) 
replacing V(ip), U(ip), G(p) and the radiation currents Jn,J v saturated by the contribution 
of the quantum mechanical mode, j%,jl- The resulting rolling force in the leading order of 
the slow roll expansion becomes the sum of the force induced by the effective action, F eG , 
and the quantum mechanical force, F q , F = F cS + F q , 

2V<*U?-ir*Vf 1 / a d f N \ 

- G c$ Uc $ + 3(^)2 > * U + 3U$[ J * v3n 2d dt)- [ } 

7. Radiation currents and their effect on inflationary 
dynamics 

The role of F eS in the inflationary evolution has been studied in ||. For the no-boundary 
and tunneling states in the one-loop approximation it equals 



X m 2 P {l + 5) (,p 2 



Fn\t = £-7^ H1T3+ O(Vier) (7-1) 



and leads to different conclusions. The no-boundary peak is realized for 1 + 5 < 0, therefore 
the point ipj is an attractor - quantum terms in F eS lock the inflaton at its constant initial 
value and give rise to infinitely long inflationary scenario with exactly DeSitter spacetime. 
In the tunneling case, the probability peak exists in the opposite range of the parameter 
( |2.4|) , 5 > — 1, and the rolling force has the quantum term which initially doubles the 
negative classical part. Therefore, the inflaton decreases under the influence of this force, 
and the tunneling state generates a finite inflation stage with the estimated e-folding number 
N ~ 487r 2 In2/A leading to the estimate on A, A < 5.5. As we shall now see, the quantum 
mechanical force F q does not qualitatively change these predictions. 



To obtain radiation currents in ( |6.5|) we expand the first order variations of ( |6.1|) to the 
second order in perturbations (A,ip,5ip) and make their quantum averaging with respect 
to the initial state. For this purpose we, first, need the initial reduced density matrix of 
the inflaton field p(p, ip') in the non-minimal model and, second, the expressions for the 
Heisenberg operators AQ(rj) in terms of quantum initial data, 5<p = S<p, p$ v = d/id(5<p). 
In |12|| it was shown that for the model with a big |£| the initial density matrix describes 



practically pure quantum state and expresses in terms of the distribution function p(ip, ip') ~ 
p 1 l 2 (ip)p 1 l 2 (ip') ) |£| ^> 1. The effective pure quantum state, ^NB,T(f) — PtvbtIv 9 ); m the 
vicinity of the probability maximum, which is located at ifi, can, thus, be approximated by 
the gaussian packet of small quantum width A - the square root of ( |2.3| ) . 

The operators of quantum disturbances in the (^-representation can be found by using 
the results of Sects. 4-5 obtained for the minimal model because it can be regarded as the 
Einstein frame for the non-minimal one. The Einstein frame for (|6.1| ) with U = m 2 p/16n and 
G = I arises by a special conformal transformation and reparameterization of the inflaton 
field |i~3| : (g^jip) —> (g^ u ,ip). We denote the objects in the Einstein frame by bars and 
identify them with those of the minimal model considered in Sects. 4-5. In this way we 
reduce all the calculations, Cauchy data setting, gauge fixing, reduction to the physical 
sector, etc. to those of the minimal model. This makes the further calculations of radiation 
currents straightforward. Below we separately consider the beginning of the inflation epoch, 
t = 0, and late stationary stage of inflation. 

At the onset of inflation t = the radiation currents equal [Tl[ 



J ^ (0) _ 384tt 2 |£| 2 [f 3 ; J' J " (0) _ 967r 2 |£| 2 [f 6 j j ' (7 ' 2) 

and (a/a)dj q N /dt{Q) = 0, where / = (A/167r 2 ) 2 |£|/|l + S\. These quantities are strongly sup- 
pressed as compared to their classical values, by a very small factor A/|£| 2 ~ AT 2 /T 2 ~ 10 -10 
related to the CMBR anisotropy |7|]. Their sign crucially depends on the magnitude of 
the parameter /, which in our model is likely to be very big, / 3> 1 (in view of the es- 
timate iV > 60 on the e-folding number and the value of |£| ~ 10 4 [0]). In this case, 
the terms proportional to / ~ 1/A 2 , generated by the kinetic terms of the radiation 
currents, (AQ'AQ'Y dominate and, in particular, lead to the positive energy density, 
£ q (0) = — Jtv(0) ~ m A p\ 2 \l + <5|/647r 2 |£| 3 <^ mp, and, as one can check [jnj, negative pres- 
sure p q (0) = — e q (0) (DeSitter equation of state). Interestingly, the sign of the quantum 
rolling force due to the homogeneous mode is independent of the magnitude of /, F q (0) ~ 
(A 2 (^f/36)//967r 2 |^| 3 > 0. In view of the expression for / the actual magnitude of this force 
is again much smaller than its classical counterpart F q (0) ~ |-P 1class (0)|/|^| <t^ |F class (0)|. 
Therefore, for the tunneling state it gives a negligible contribution to the full rolling force. 
For the no-boundary state, the initial effective force (|7.1|) vanishes, but the only effect that 
the positive F q (0) can produce is that it shifts the equilibrium point from ipj to slightly 
higher value, at which again the system undergoes endless inflation. 

At late stationary stage of inflation the dynamics of the classical background can be 
approximated by the ansatz a = coshfif (ip)t]/H (ip), if ~ ipj, and the resulting radiation 
currents are dominated by the contribution of the growing mode q + ~ smb. Ht/3a(t)(f) — > 



const, Ht 3> 1: 

JiV " 4 8647r 2 |e| 2/ ' ^" 2 8647r 2 |er' l ' 

Similarly to the onset of inflation, they are strongly suppressed relative to the classical 
values by the factor A/|£| 2 ~ 10~ 10 . The energy density of the quantum mechanical mode, 
e q = —j q N ~ — A 2 1 1 + <5|mf>/54(167r 2 ) 2 |£| 3 ), \e q \ <C rrip, is negative, while the pressure 
is positive p q = —e q |lTJ (anti-DeSitter case). Apparently, this is a manifestation of the 



ghost nature of the mode q whose kinetic term enters the action (|5.1| ) with the wrong sign. 
Radiation currents (|7.3|) lead to the quantum rolling force ( |6.4[ ), F q ~ 0, which vanishes in 
the leading order of the slow roll expansion. 

8. Conclusion 

The dynamical contribution of the quantum mechanical mode to effective equations 
turned out to be very small - it is strongly dominated by the effective rolling force. This 
property was actually conjectured in ||, and now it is quantitatively confirmed. Thus, 
this mode cannot change the dynamical predictions in spatially closed model with strong 
non-minimal coupling. As a model of the low-energy quantum origin of the Universe only 
the tunneling state remains observationally justified, because the no-boundary wavefunction 
generates infinitely long inflationary stage. The role of this mode should not, however, be 
underestimated, because its effect is model dependent, and might be important in other 
models generating initial conditions for inflation. Moreover, the inflaton mode simulates the 
DeSitter and Anti-DeSitter effective equations of state, e+p = 0, respectively at the onset of 
inflation and at late times. The sign of its energy density contribution can change depending 
on the balance of the potential and kinetic terms of this ghost mode. Therefore, it is not 
quite clear at the moment, what can the role of this mode be at post inflationary epoch. 
A natural question arises if this mode can be responsible for the present day observable 
acceleration of the Universe as an alternative to quintessence or be capable of inducing 
DeSitter-anti-DeSitter phase transitions in cosmology? This question is subject to further 
studies. 
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